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INTRODUCTION 

In this meeting we have talked mainly about nuclei and baryons but I will concentrate on 
mesons. The reason is that the lightest mesons are the simplest bound states present in 
Quantum Chrome Dynamics (QCD). They have the minimal number of constituents and 
are the lightest state so should be spatially the simplest. Their properties are to a large 
extent determined by Chiral Symmetry which enforces vanishing masses in the limit of 
zero current quark masses, the chiral limit, as well as vanishing interactions in the zero 
momentum and chiral limit. It is the combination of these two properties that allows us 
to produce a well defined low energy theory, Chiral Perturbation Theory (ChPT), as a 
consistent approximation to QCD. In the remainder I review ChPT at two loops. 



(Effective) Field Theory, ChPT and lattice QCD 

The underlying idea is really the essence of (most of) physics. Use the right degrees 
of freedom. When there is a gap in the spectrum with a consequent separation of scales, 
we can build the theory containing only the lighter degrees of freedom and include 
the effects of the neglected high mass/energy states perturbatively by building the most 
general (local) Lagrangian with the low mass degrees of freedom. This leads in general 
to an infinite number of parameters so no predictivity is left. But when these terms can 
be ordered in importance by some principle, usually referred to as power counting, we 
have a finite number of parameters at any given order and thus an effective theory. 

We need to use field theory since it is the only known way to locally combine quantum 
mechanics and special relativity. A Taylor expansion cannot be used because in the 
chiral limit the continuum of states gives it zero convergence radius. As important is 
the fact that off-shell effects are fully under control. The freedom allowed by this is 
fully described by extra free parameters. In addition it is systematic, all effects at a 
given order in the expansion are included and errors can be estimated. Drawbacks are of 
course the large number of free parameters, but do not forget the fact that while a model 
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Mesons Quark Flow Quark Flow 

Valence Sea 




FIGURE 1. An example of a meson (thick line) loop diagram (left) with the quark diagrams (right), 
there are gluons all over in the right diagrams. The Valence (arrow line) and Sea (dash arrow line) quarks 
can be treated independently in lattice QCD. 

can have few parameters there is a large freedom in the space of possible models hidden 
behind it, and as always the expansion itself might also not converge. 

For ChPT the included degrees of freedom are the Goldstone Bosons from spon- 
taneous breakdown of chiral symmetry, identified with the pseudoscalar octet, n"^'^, 
K^,K^, and r\. The powercounting principle is dimensional counting, counting pow- 
ers of generic momenta p, and the expected breakdown scale is the mass of the neglected 
resonances, on the order of Mp. 

Chiral symmetry is the interchange of m, d, s quarks which in the limit — > can be 
done independently for the left and right chirality since the QCD Lagrangian, 

■^QCD = ^ [iqtlpgL + iqRlpgR - mq {quqi + qiqu)] , ( 1 ) 

q=u,d,s 

does not couple left and right in the chiral limit. The symmetry group is SU x 
SU{3)r and is spontaneously broken to the diagonal subgroup by {qq) = {qLqR + 
qRqi) 7^ 0. The 8 broken generators correspond to the pseudoscalart octet. ChPT in 
its present form was introduced by Weinberg, Gasser and Leutwyler [1, 2, 3] and 
introductory lectures can be found in [4]. 

In the case of lattice QCD we need to extend this. One distinguishes between valence 
and sea quarks and they can be treated independently. A simple example of how a meson 
loop contains quark loops is shown in Fig. 1. In general it is much easier to get to light 
valence quarks than to light sea quarks. A more extended discussion can be found in 
[5, 6]. The version of ChPT with sea and valence quarks separately is called Partially 
Quenched ChPT (PQChPT). 



TWO LOOP: GENERAL 

The Lagrangians needed at the first three orders, p^, and p^, are known and the 
number of parameters and their notation are shown in Tab. 1 . The parameters at order p^ 
go back very long time. The classification and infinities at p^ can be found in [2, 3] and 
at p^ in [7, 8]. The knowledge of the infinities provides a very useful check on explicit 
calculations. 
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TABLE 1. The number of free parameters and their names for the 
first three orders in mesonic ChPT. The numbers refer to the number of 
low energy constants and unmeasurable (high energy) constants. 







2 flavor 




3 flavor 


3+3 PQChPT 




F,B 


2 




2 


Fo,Bq 2 






7+3 




10+2 


L\,Hl 11+2 


p' 




53+4 




90+4 


Kr 112+3 



The replica method allows to obtain PQChPT from the np flavor case of [7, 8]. 
Important is that the 3 flavor case is a limit of the 3+3 PQChPT one. Determining the 
parameters in PQChPT allows to obtain those of 3 flavor ChPT via the Cayley-Hamilton 
relations of [7]. 



TWO LOOP: 2 FLAVOURS 

In this case most calculations have already been performed. 

• rr^;rV[9]. •Fv^{t),Fs::[\2l 

• yy ^ K~ , Ftj:, nij^ [\Qi\. mn^lvy[\3\. 

• ;r;r-scattering, Fn,mTi [11]. 

Here there is in general a rather good convergence and it turns out that for many thresh- 
old quantities the values of the cj. are not numerically important. These calculations are 
now often combined with dispersive methods for very precise theoretical predictions. 
An example is the full description of TtTt scattering [14]. 



TWO LOOP: PQCHPT 

This subject is only beginning. The charged pion mass for the case of all valence 

masses equal and all 3 sea quark masses equal is published [15] and the result for the 
decay constant in this case will be published soon. Planned future work includes all the 
necessary mass combinations as well as the two sea quark case. 

In the actual calculations heavy use is made of the symbolic manipulation program 
FORM [16]. The major problem is simply the sheer size of the expressions involved in 
the more general mass case due to the appearance in PQChPT of many partial fractions 
of differences of quark masses. An example of results is shown in Fig. 2. 



TWO LOOP: 3 FLAVOR OVERVIEW AND nn, nK SCATTERING 

Many calculations have been performed here. A (to my knowledge) complete list is 
given below, where the quantities which are determined from this calculation are shown 
in brackets. The symbol 11/,/^ means a two-point Green function of currents / and / with 
the quantum numbers of meson M. 
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FIGURE 2. First results from PQChPT at two-loop level [15]. Shown is the relative correction to the 
lowest order mass squared at order p'^ (left) and + (right). The valence quark mass is given as 
Xi ~ IBQiUgv and the sea quark mass as X4 = ^Bom^s with X4 — Xi tan 9 (curves labeled with angle 9) or 
X4 — 0.125 GeV^ (A). Xi corresponds to the lowest order meson mass squared. The boundary for ChPT is 
reached for some mass squared around 0.3 GeV^ 



• Fn,K,-n, m^^K,n {ruu ^ nid) [23] {U^ -j^^.tnu/nid) 
We now perform a general fit to experiment and check how well the whole system works 
and determine as many parameters as possible. First we need to identify a series of 
basic inputs to determine most of the parameters. The procedure is described in detail in 
[22, 23]. The actual inputs used are 

• K^^. F{0), G(0), A from the E865 BNL experiment [31]. 

• m^Q, m^, m^+, m^g Electromagnetic corrections including quark mass effects[32]. 



• The C[ as estimated using single resonance approximation. 

These fits are performed varying the resonance input used for the Cf by varying the size 
by an overall factor of two and the scale at which the saturation is applied. The final 
fitting result with errors is shown for some typical inputs in Tab. 2. The C[ can be de- 
termined from experiment in several cases. When the dependence is purely kinematical, 
e.g. curvature of a form-factor, it works reasonably well [12, 25]. For those with a mixed 
quark mass-kinematical dependence, e.g. the slope of f+{t) in Kfi, it works OK [26]. 



• UvvK,^vvn [17, 18] 

• UvvK [18, 19] 

• ^AAn, HaAt], Fj^,ri, tn,^,r], [18, 20] 

• nss [21] (L^,L^) 

• UvVK, ^AAK, Fk, niK [18] 
.i^M[22](L'i,L^,L^) 



•Fvn,FyK+,Fyj,o [24, 25] m) 

• Ka [24, 26] (Vus) 
•Fsn,FsK[21] {L\,Ll) 

• KK [29] 

• kK [30] 




92.4MeV,Fj^+/F; 



1.22 ±0.01, 2m,/(m„ + mrf) =24 or 26. 
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TABLE 2. The results of the fits for several input values of ^-4,^5, see text. 





fit 10 


4 

same p 


fitB 


fitD 





0.43 ±0.12 


0.38 


0.44 


0.44 




0.73 ±0.12 


1.59 


0.60 


0.69 


lO^L^ 


-2.53 ±0.37 


—2.91 


-2.31 


-2.33 




= u 


= u 


= U.J 


= u.z 




0.97±0.11 


1.46 


0.82 


0.88 


10% 


= 


= 


= 0.1 


= 




-0.31 ±0.14 


-0.49 


-0.26 


-0.28 


10% 


0.60 ±0.18 


1.00 


0.50 


0.54 




0.736 


0.991 


1.129 


0.958 


2 4 6 


0.006,0.258 


0.009,= 


-0.138,0.009 


-0.091,0.133 


2 4 6 


0.007,0.306 


0.075,= 


-0.149,0.094 


-0.096,0.201 


2 4 6 


-0.052,0.318 


0.013,= 


-0.197,0.073 


-0.151,0.197 


mu/nid 


0.45±0.05 


0.52 


0.52 


0.50 


Fo [MeV] 


87.7 


81.1 


70.4 


80.4 


^k- „4 6 
F^-P 'P 


0.169,0.051 


0.22,= 


0.153,0.067 


0.159,0.061 



The remaining ones are difficult to estimate, the question here is what type of scalar to 
use, we know it is not the sigma [12, 27] but what else is an open question. 

The whole procedure is repeated for a range of input values of ^5. Some examples 
of the fits are given in Tab. 2 where we quote the main fit (labeled fit 10), the same one 
but at order rather than as well as two with L\,Uf^^ 0. Fit B is one where all scalar 
form factors behave nicely [27] while fit D is the one where the threshold parameters for 
KK and kK scattering are well fitted as discussed in more detail in [30]. Fig. 3 shows the 
constraints from TtK and kK scattering as well as the region of fit D. 



TWO LOOP: Kb AND Vus 

The CKM matrix element Vus is a fundamental parameter we like to determine as 
precisely as possible. In addition, the relation 1^^^+ iKvP = 1 is broken at about 
the two sigma level by the experimental values quoted in the PDG 2002[33]. It is thus 
important to update both theory and experiment. The theory [34] has been updated in 
several ways. The photonic corrections have been properly calculated in the modern 
ChPT language [35] and the form factors f+{t), /- (t) are known to p^ in ChPT [24, 26]. 

The experiments have been mainly analyzed by assuming a linear parameterization of 
the form-factor. This is not sufficient as was pointed out in [26]. The value of \f+{0)Vus\ 
changes by 0.9% (0.6%), data from [36] ([37]), when a linear fit is used as compared 
with a quadratic one where the curvature is determined from ChPT [26]. The newer 
experiments [38, 39] have now detected the curvature and find values in agreement 
with the ChPT prediction of [26]. The measured branching ratio in both the neutral 
and charged channel has also increased, [38, 40, 41]. Both combined lead to an increase 
in the value of |/-|-(0)V„j| and using the value of /+(0) from [34] the unitarity problem 
is resolved. 
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Ttjt constraints 



jtK constraints 




4 -0.2 0.2 0.4 0.6 



FIGURE 3. The constraints on L\,Ll from KK (left) and kK (right) scattering [29, 30]. 



The remaining problem is to accurately predict the value of /+ (0) . In [34] one-loop 
ChPT and a quark model estimate of the higher orders were used. The calculation 
gives loop corrections of about one % to the one loop result [24, 26]. The full analysis 
including updated values of all the inputs is in [26] but a more important point was made 
in [26] as well. The low energy constants that appear in the value of /+(0) can be 
determined experimentally from the scalar form factor /o(0 since [26] 



hit) 



8 



8 



m 



m. 



iFK/F^-l)-—^t'C\^ + ^{t)+m■ 

^7T 



(2) 



In this equation everything is known except the values of €[2 and C34, and correlations 
between Fjc/F^ and V'tu[42]. Since /+(0) = /o(0) a measurement of the slope and 
curvature of the scalar form factor allows to determine /+(0). Dispersion theory allows 
to relate some of these quantities as well. A first analysis [43] leads to an estimate that 
essentially cancels the pure loop contribution yielding a very small total p^ correction. 



CONCLUSIONS 

The two flavor case in ChPT at two loops is an almost finished subject. The three 
flavor case is in progress. Many calculations have been done and things seem to work 
but convergence is sometimes slow. L4,Lg are nonzero but reasonable for large Nc 
expectations. An example of clean predictions even with the many p^ constants is the 
case of Kit, decays and the determination of Vus. The partially quenched case, relevant 
for lattice calculations, is just at its beginnings at two loop order. 
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